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I. INTFODUCTIO:^ 

The naxlmuii power level of a nuclear reactor is fre- 
quently limited by the rate of heat removal. The actual de- 
sign of a power reactor represents a coit^sromise between the 
requirements of heat transfer and those of nuclear physics. 
Nuclear reactor theory is not sufficiently accurate to guar- 
antee that a particular reactor (especially one of a new 
design) will go critical. Therefore, the nuclear characteris- 
tics of a reactor design must be obtained before construction 
of the full scale reactor, Sono of this information can be 
obtained from the us© of a subcritlcal assembly « 

Experiments with uranium graphite subcritlcal assemblies 
can be carried out to determine the various nuclear constants 
of a proposed reactor lattice. These lattice constants are 
the material buckling, multiplication constant, lattice dif- 
fusion length, and thermal utilization of the imlt cell in- 
cluding the process tube assembly. 

In this investigation six lattice configurations were 
ccnnldered, and the buckling of each was determined experi- 
mentally, lattice colls of 6 in., 0,5 in., and 12 in., de- 
noted I, II, and III respectively, constituted the throe geo- 
metric lattice arrangements. For each of these, both the 
”wat” and *'dry” cases woi'e considered. The measurements and 
calculations with i-iater present in the cooling annulus were 
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designated as "wot”, and those with the process tube, but 
without water, as ”dry”. 

These six lattice configiorations were also considered 
in the theoretical analysis. The buckling of each of the 
six lattice oonf igurat ions was determined theoretically by 
two methods. Only one size f\iel element and one size proc- 
ess tubing was used in the experimental and theoretical 
analyses of this investigation. 

The methods of Murray (1) and of Rumsey and Volkoff (2) 
(3) were used to determine theoretically the thermal utlllza 
tion and the resonance escape probability. Hurray's method 
was shown to give loss conservative results than the method 
of Fuiasey and Volkoff. The theoretical multiplication con- 
stant derived for each conf Igiiration considered was based 
on literature and calculated values of the various diffusion 
lengths, the Fermi age, and the various raaox’oscopic cross 
sections. The product ot and was assumed to be es- 
sentially constant for the various configurations investi- 
gated. 

Analyses of the data for the given configurations were 
carried out by use of the classical formula for the buckling 
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A cor^arlaon between the **wet” and ”dry‘* cases of the 
throe lattice cells considered was presented using the experi- 
Biental and theoretical values of the buckling* By an analysis 
of the results obtained in this investigation, tho effect of 
the water coolant on the nuclear properties of the various 
lattice configurations was detortnlned* 
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II. REVIEW 0? LITERATURE 

The theoretical calculation of the multiplication factor 
for a heterogeneous reactor with coolant la an extension to 
the method of calculating the thermal utilization and reso- 
nance escape probability factors for an infinite heterogeneous 
uranium-graphite reactor presented by Weinberg ( 4 ) » i$) • 
Guggenheim and Pryoe (6) extended the method to include a non- 
absorbing coolant and a thin sheath of absorbent aluminum 
around the fuel element. Houston ( 7 ) extended the method to 
include neutron absorption In a coolant. Houston's results 
arc more suitable for gas or heavy liquid notal cooled re- 
actors due to his assumptions concerning the coolant 's negli- 
gible moderating properties and small volume fraction of tho 
lattice coll. Rumsey and Volkoff (2), ( 3 ) extended the above 
methods further to inclxxde tho moderating effect of a water 
annulus . 

ffurray (1) presented a simplified method of calculating 
the thermal utilization which considers all but tho tiranium 
and the graphite in the lattice coll as poisons which reduce 
the basic uranium graphite thermal utilization. Kui’ray's 
method for the calculation of the resonance escape probability 
does not consider the so called poisons. The Rumsey and Vol- 
koff method appears to bo at the present time the most com- 
prehensive theoretical analysis of a water coolod uranlxim- 
graphlto lattice coll. 
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Th« experiments with aubcritic&l aseeabliea of th© Han- 
ford Atomic Producto Operation included an effort to determine 
the effect of a water annulus around the fuel element on th© 
buckling, The effect of alualnum, air, end other impurities 
in the lattice upon the buckling was also inveotlgated (8), 

A few Eieasurcments were taken using internally cooled cylin- 
drical fuel rods, but the laajority of the coolant oinulation 
experiments were with solid cylindrical fuel rods of natural 
uranium surrounded by an annulus of vrater# The objectives of 
th© subcrltlcal assembly exp©rl?«cnts at Hanford (3) were to 
detorralne oxperlrcontally the various lattice constants. 'Fhe 
existing lattice theory was supplemented and improved whore 
possible. Irfittlce measurenonts were taken with three differ- 
ent fuel rod diameters: 0.925 in., 1.175 in., and 1.36 in. 

The lattice spacing and th© water- aluminum-uranium ratios 
were varied for each slug siae. 

The results of the series of tests using both "wet" and 
"dry" lattices vfere ccmpllod by Clayton (9) . A cojigpariaon of 
theory and experimental roaults for one fxiol element diameter, 
six lattice spaclngs, and the "wet" and "dry" configurations 
for the Hanford Atomic Products (^©ration was done by Cast 
and others (2) . Cast used th© ir^thod of Humsoy and Volkoff to 
calciilat© th© theoretical values of the thermal utilisation 
and th© resonance escape probability. 

Related conclusions (8) from the Hanford work with water 
cooled uranium graphite lattices were: 
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(a) Por a givea natural uranluw fuel assembly, 
there is a mxlmum in the buckling vs. C/U 
atom ratio curve. The slug diameter and the 
water annulus thidcncss arc inversely pro- 
portional to the value of C/U at which the 
maximazB occurs* Tbia maximum occurs between 
C/U values of 50 and 100. 

(b) The buckling “dry” is greater than the buck- 
ling "wet" for larger lattice spaolngs (g/D 
atom ratio). As the lattice spacing Is de- 
creased, the buckling curves cross and the 
”wet” buckling becomes greater than the "dry” 
buckling. For this lower c/U atom ratio 
region, externally cooled slugs exhibit a 
"fall-safe” behavior on the loss of water. 

The following conclusions were a result of the Hanford 
experiments(2) . The state of the hetorogenoous graphite 
tiranlum lattice theory at present is not completely satis- 
factory. More experimental weanuroments of individual param- 
eters such as f and p are necessary in order to establish the 
basis for an improved theoretical analysis. Corrections for 
variations In reactor material densities and purities, and 
allowances for holes, such as for control rods, need further 
Investigation (2). Finally, the treatment of the energy spec- 
trtjm of neutrons, fast to thermal, reeulros more accurate 
development of neutron resonance and capture reactions. 
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III. TB?:0HE7ICAI. ANALiaiS 
A. General Procedxire 

The product of the fast fission factor, 6 , the primry 
fission neutron factor, »| , the thermal utilisation factor, 
f, and the resonance escape factor, p, is the Infinite nulti- 
plication constant, kj, « 



€Vjfp SB (2) 

A theoretical analysis of each of the factors would have 
been extremely conplicated without certain simplifying as- 
siDiptiona. Theiie aoauiaptiona arc stated where applicable in 
the following theory. Each of the factors was considered 
separately. The product 6V| >«is determined and assumed con- 
stant for all the lattice conf Iguratlons considered. Expres- 
sions for f and p were determined using the methods of Murray, 
denoted Method A, and the methods of Humscy and Volkoff , de- 
noted Method B, for the ”wet" and ’’dry” cases of Lattices I, 

II and III. 

Using values from the literature for the lattice diffu- 
sion length and the Fermi age with the theoretically deter- 
mined infinite multiplication constant, the buckling was deter- 
mined from Fqmtlon 1» The theorotlcal buckling was determined 
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for Lattices I, II, and III for both the ”wet" and "dry” 
cases using both Methods A and B. 

The affect of the air annulus and air holes on the neu- 
tron densltj' in all lattice geometries was neglected due to 
the magnitude of the effect in comparison with other material 
effects in the lattice. 

Subscripts used in this analysis were defined as follows: 
u natural uraniuis 
al dug can, 2S aluminum 
w water 
g graphite 
m moderator 

p process tube, 6lS aluminum 
a absorption 
8 scattering 
re resonance 

aa scattering "absorption” 
f slow neutron fission 
t transport theory 

Any deviation of notation from the above list was individual- 
ly defined. 



B. Past Fission F-’actor 

The fast fission factor, € , is the ratio of the total 
number of neutrons produced by fission to the number of neutrons 
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produced by flssloa by therrnal neutrons. Assuming that the 
primary neutron source is distributed approximately uniformly 
over the fuel, a slr^lifled version of the equation given by 
Claestone and Edlund (5* p* 278) for the fast fission, factor 
is 



£ - 1 



(j/ - 1) - 



1 : 



0~c 1 ^f 



1 - P 



(- 



j/CTf 
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where P is the probability that a fission neutron born in 

a rod will make a collision inside the rod in which 
it was created} 

the absorption cross section for fast neutrons, 
that is, the sum of the fission and the capture cross 



sections} 

0". is the elastic cross section for fast neutrons} 
is the Inelastic cross section for fast neutrons; 
and CTt is the total cross section for fast neutrons. 

Using fast neutron cross section values as given by Gug- 
genheim and Pryce (6, p. 5l)» Equation 3 reduces to 



0.09liO P 
1 - 0,521 P 



€ - 1 => 



(4) 
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These nuncrlcal coefficients are such that when the fuel elo- 

laont raditis is loss than 1*7 cirj the calculated value of ^ 

agrees reasonably well with oxperlmontal results obtained by 

Argonne National Laboratory (6, p*5>D • For natural xiranlusa 

of density 19 g/cm"^ and a rod radius of 0,5 in., the value of 

Her, r. is 0.262, where M is the nuiaber of uranium atoms per 
t u 

unit volume, and r is the radius of the fuel rod. The quari- 
tity P is a function of IfCT. r . From a relationship of P and 
N OTj. r^ shown by Guggenheim and Pryce (6, p, 51 ) and the above 
computed value of N O"^ r^, the value of P determined vms 0.25* 
Solving rquation I4. for € gave a value of 1.027. This 
value for the fast fission factor was assuned constant for all 
conf igurationa investigated. 

C. Primary Fission Neutrons 



The number of fission neutrons released per thermal neu- 
tron captured in the fuel la the factor, kj . Glaastone and 
Edlund define Kj (5» p. 83) as 







(S) 



where ^ is the average number of neutrons produced per thermal 



neutron fission; 
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Is the total cross section for absorption of thermal 
neutrons. 

The microscopic thermal fission crosa section for natural 
uraniuia is barns. Applying the maxwellian distribution 

correction of 1*128 and the "not - 1/v” factor of 0.981 for 

to this value gives an effective fission cross section of 
3*635 barns for natural uranium. The microscopic thermal ab- 
sorption cross section for natural uranium is 7*68 barns. Ap- 
plying the maxwellian distribution correction and a "not - 1/v” 
factor of 0.990 for natural \iranium to this value determines 
an effective absorption cross section of 6.74 barns for natu- 
ral uranium (1, p, 32). 

Murray (1, p, 98) gives a value of 2.46 for ^ • Sub- 
stituting the above quantities in Equation $, the calculated 
value of trj was 1,327. 

It is normal practice in working with subcrltical assem- 
blies to experimentally determine k», and then divide by 
the product t pf to arrive at a value of . The value of 
found by subcrltical assembly experiments was 1.308 neutrons 
per thermal neutron captured (10, p. 65)* This value la in 
reasonable agreement to the computed value of 1,327. The val- 
ue of irj used in the following analysis was 1,303 and was as- 
sumed constant for all configurations investigated. 




M gl| III I «4 



«^4§*^ 1 % 2 . Im / c» 

^ 4yb4i4«^^« -ii<4te »«|««H» 4i «-t«««i^ |^9an«B» ^ •* 




^•V l^>C^ •' ««i •* «BSl •> •** « 



• % ami 1 1 1 1 tf» ( 



j 



12 



D, Thoi'iaal Ut ill sat ion 

In heterogonoous lattice theory the finite lattice re- 
actor is assumed to bo an infinite lattice array. Tho prob- 
lem of calculating* f and p Is sinpllfiod bj replacing tho 
square colls by equivalent cylindrical unit cells. Table 1 
gives the equivalent cell dimensions for Lattices I, II, and 
III. Assuming that little ebango in the overall neutron 
flux is experlencod in traversing on© cell, an individual 
coll nay bo roprosontative of tho total. 

The methods of Murray {!) and of Ktimeey and Volkoff (2), 
(3) wore used to determin© theoretically tho thcxnnal utiliza- 
tion. 'Ihese niethods were considered soperatoly and were de- 
noted Methods A and B respectively. 

1. Method A 

The depression of the thermal neutron flux in tho fuel 
and adjacent tnodorator complicates the calculation of the 
thei'mal utilisation. Average flux and volxirie weighting fac- 
tors must both be Included in the fractional absorption of 
the fuol. Prom tho spatial vai'lation of the flux in a coll, 
the average values of the fluuc cen be determined. 

The thermal neutron diffusion equations for the la'aniiim 
and the moderator from diffusion theory are 



D sf - 0 Z+s-0 

U u u 



(6) 
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Equations 6 and 7 can b© sl’^tpllfled «lth the following aa- 
•unptiona: 

(a) The variation of neutron flux at a given 
coll radius will be the same everywhere, 

(b) There are no thermal neutrons produced 

in the uranium, is zero. 

(c) The production rate of thermal neutrons 
in the moderator is independent of 
position, 3 is constant. 



(d) The flux is assumed not to vary along 
the fuel coll axis. 

Using the above assumptions and cylindrical coordinates. 
Equation 6 for the fuel becomes 



vdier© X » ^r and is the inverse diffusion length of 
uranium, and Equation 7 for the moderator becomes 




(8) 
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where P la the diffusion coefficient for the moderator. 

4f la the aquere of the inverse diffusion length. 

TtX 

Equation 8 for the fuel is the modified Bessel equation 
of zero order. Equation 9 for the moderator la the inhomo- 
geneous type of the modified Bessel equation of zero order* 
The solutions of Equations 6 and 9 are 



= A Iq iifyT) 



( 10 ) 
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respectively, 7ho boundary conditions used in the evaluation 
of the constants A, C, and C are neutron flux and current 
continuity at r = r^ (urariium rod radius) and the condition 
that » 0 at r a Tg (the equivalent outer radius of the 

graphite coll). 

Applying the condition of zero neutron current at r^ in 
Equation 11 
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thus the moderator flux. Equation 11, becomes 
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Continuity of flux and current ct r^ give* a pair of ©qua- 
tlone from Equations 10 and 13 which can be solved simul- 
tanooualy to give 



A ^ CT ml ra u 



(lu) 



and 



S_ 



\>^uh 



(IS) 



)a 



whore 



«1 .K, (rf,r) - (dP^r) 

X m 2 



and 






& s 


I 

ra 


XA m ' 


with 




dH^Cx) 






dx 




The 


thermal 



”a ''a Io< ^u*-u> «1< '<'a'-u> " r>u '''u ^u-u> «0< ^a*-u> 



- Hj^(x) 



neutrons absorbed in the fuel to tho total thermal neutrons 
absorbed* This ratio is equivalent to tho ratio of the thor- 
mal neutrons absorbed in the fuel to the total thermal neutron 
aovu’c© assuming no net leakage from the cell* 
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Gocfcinlng T'^quation 10 and 16 Into 



f =» 



ff) ^ /“ ^0 

'a' ra 0 
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and solving the integral ( 1 , p« 305 ) gives 
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Corabining Equations II 4 . and I 8 and solving Tor the reciprocal 
of f gives 
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= 1 + 
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P + (E - 1) 
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where 
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P 1» the dlsadvanteige factor of the uraniua wiiich Is the 
ratio of the neutron flux at tho rod surface to the average 
fltix in the rod interior (S>» p«270). 



p » 




(22) 



Tho average thernml neutron flux in the moderator Is greater 
than the thermal neutron flux at the rod surface because the 
diffusion coefficient in the moderator is finite. The ad- 
ditional absorption is measured bj the quantity (E - 1), the 
excess absorption term. 

Murray (1, p. lOD) assumed that any poisons that are 
tolerable do not appreciably disturb the basic noutron flux 
distribution from that of the configuration with fuel and 
moderator only. Therefore, the aluminum cladding, coolant, 
and tho process tube assembly can all be treated as poisons. 
The reciprocal thermal utlllaation with this isodlf ication can 
be written 




n poisons 

h A 

L ^ K 



1*0 



(23) 



where the average neutron flux in the added poison or absorber 
is Fi- 



1*9 
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Assuming that for the cladding, water annulus, and proc- 
ess tube the ratio 0 ^/ is apiiroximatoly equal to 

solved directly. Ilurray (1, 

p. 100) suggested that a slight refinement to Eouatlon 23 may 
be made by letting be the difference between the actual 
poison cross section and the moderator cross section that it 
physically replaces. 

Lit tier and Raffle (10, p. 94) give a useful approxima- 
tion for the excess absorption term, E - 1. 



E - 1 * 
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The constants used In Equation 23 were listed in Tables 1 
and 2. The equivalent cell dimensions were derived from di- 
rect Hsoasuremeats end calculations. The thermal neutron macro- 
scopic cross sections for aliiminum (cladding), water, and 
graphite were calculated from Equation 25 using literature 

values of cr (1), (5), and (11). 
a 

I = -£f^ (25) 

where *?q la 6.02 x 10^^ nuclel/mol©, 
and A is the mass number. 
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for the eliinlnum process tube was calculated using a 
modification of Equation 25 to account for the alloying ma- 
terials in the aluminum alloy. The laaxv/elllan diatribatlon 
correction and a ”not-l/v” correction was used to determine 
an effective absorption cross section for natural uraniiua. 

All the macroscopic cross section values used in this theo- 
retical analysis are shoim in Table 2. 

Du© to the limitations of simple diffusion theory, trans- 
port theory must be used for the calculation of J!^^. Kurray 
(1, p. 88) gives for the uranium fuel the relationship 




tanh 
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( 26 ) 



^dioro is the total cross section. 

Solving Equation 25 using cross section values from Ta- 
ble 2 gave a value of 0.675 cm“^ for The inverse dif- 

fusion lengths for graphite and water were found using Equa- 
tion 27 and values of the diffusion lengths from the literature 

(5). 



s I * ( Z^/n)’-/^ (27) 

The thermal utilizstlons for the dry lattice configura- 
tions were calculated by omitting the term for water in Equa- 
tion 23. The volume ratios in Equation 19 and 23 are equiva- 
lent to area ratios. 
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The conotanta listed in Tables 1 end 2 were need in 
Eqiiations 19, 20, 23, end 21}. to evaluate the coariposite ther- 
mal utilizations for the various lattice configurations con- 
sidered. A Bessel function table (12) was used to evaluate 
the disadvantage factor, P, in Equation 20* Equation 2l4. was 
used to evaluate the excess absorption term, (E - 1). 

The resulting thermal utilizations coraputed by Method A 
are listed in Table 3 fof the various lattice conrigvu*ationa. 
2. Method B 

A second method of calculating the thermal utilization 
of uranium is the met od of Rurasey and Volkoff (2), (3)* This 
method is an extension of prior work in heterogeneous lattice 
theory (I|.), (6) and Includes the moderating effect of a water 
annulus in the lattice cell. 
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The "relative abBorptlon" term, denotea per therpial 
neutron cap toured In the uranluTa, the number of thermal neu- 
trons which would be captured in the jth nodiua if the neutron 
density in the jth radium were uniformly equal to the neutron 
density at the uranlxia-aluialnum interface. 
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where F is the disadvantage factor of the uranium. Equation 

20 . 

The "self blocking” terra, S^, denotes per thermal neutron 
absorbed in the uranium, the excess number of thermal neutrons 
captured in the jth raediura due to the excess density of neu- 
trons in the jth medium over the neutron density at the i - jth 
interface. 
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( 33 ) 



where (E - 1) is the disadvantage factor of the moderator or 
excess absorption term given by Equation 24» 

The "blocking tern”, denotes per thersal neutron ab- 

sorbed in uraniiua, the excess number of thermal neutrons 
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c»ptur«d in th« Jth isediua du« to neutron density rise across 
the ith nedius}. 
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(35) 



In the above six equations is the inverse diffusion 
length of water, and t„ Is the thickness of the water annulus. 
The crose section values ar© listed in Table 2, and the voluras 
values per slvig are listed in Table 1. The value of ( -1} 

is given by Equation 2I4. idiich is an approxiciation of the coss- 
plox Bessel function in Equation 21. Tlis approxinatlon is 
considered to be stiff iciently accurate for t!ie two lattice 
geometries. 

Tho production rate of tberxaal neutrons per unit volu«e 
per second for water and graphite arc denoted and c re- 
spectivoly. The ratio of the slowing down power of water to 
graphite, »qualc 20 (3, p. 22). 

W 

Piguro 1 and 7igur© 2 show tho neutron density distribu- 
tion in the dry and wet lattice conflgtu'atlons respectively. 
The several tems of the cor^etltive absorption. Equation 29, 
are shown scheiaatlcally in these diagrams. The placement of 
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tho individual absorption terms, defined in Equations 31 
through 35* and the relative macnitudo of these various terms 
can also be seen in the diagram. Any absorption of neutrons 
by the air In the air annulus was neglected. 'Vhe net flow of 
thcrnal neutrons across the boundary of a cell was assumed to 
be zero. 

Tho term S^. aafl «hloh ar. 

shown in Figure 1 and Figure 2 were asaumod as negligible and 
were not included in Equation 29* At the uranims-a luminura 
interface, continuity of the neutron flux density and current 
density was assusKsd. The neutron current density was also as- 
sumed as linear through the thin mediums (al, p, or w) con- 
sidered. 

As an example, in deriving the term, the aluminxiia 
surface neutron flux is 



V surface 
^al 




and with Equation (10) is 




By definition the ”self blocking” tern for aluminum is 
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(36) 



and taking from (10, p* 89) 



X.®''® = :r- 1, ( 



u 



^ r 1 u u 
u 



and substituting gives 



^ 1 4. 2^2 
al 2 al al 



(37) 



Similarly 






1 fc 2 jf 2 

2 w w 



(32) 



and 



S i 

P 2 p 



2 <jf 2 

P 



A calculation of Equation 37, showed that with = 

O.OO 2 I 1 .I 4 . (10, p. 93) and « 0*102 cn, that 

_ 

» 1.27 Jt 10 This value is of a nagnitudo whloh can be 
considered negligible. 

o 

The very small factor, ( t^)^» for the cladding and 

process tube c&vises the terms 8 ^]^, S^, ^^xw* ®alp* ^alg* 

B to become negligible with respect to the other terms in 
PS 

Plgiire 1 and Figure 2. With those eliminations. Figure 1 and 



Figure 2 agree with Equation 29 
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The constants listed In Tables 1 and 2 were used in 
Equations 28 tlmourh 35 to evaluate the wet conf ipiucation 
thermal utilizations by the Ruasey and Volkoff method. Eo.ua- 
tion 20 was used to evaluate ? In Equation 31* The value 
(E - 1) in Eqimtlon 33 was evaluated by using Equation 24« 

The thermal utilizations for the dry conf Ig’orations wero cal- 
culated using Equations 28 through 35 with the terns for wa- 
ter omlttad. The resulting thermal utilizations, hothod E, 
are listed In Table 3 for the various configurations. 

E. Heeonance Escape Probability 

The fraction of fast fission neutrons that reach thermal 
energy witho^it ©xporienclng capture while slowing down, is 
oalled the resonance escape probability, p. The methods of 
rtoray (1) and of Ruiasey and Volkoff (2) were also used to 
theoretically doternine the resonance escape factor. These 
aothoda are considered separately and are denoted Methods A 
and E respectively. 

1. Method A 

Hurray (1, p. 92) gives the hetero„encous equivalent for 

p as 

\ 4 If r.” 



p « exp 



(38) 
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whor© is the effcotlvo i^croocooic cross section. 

I'e 

The form of Equation 38 nu£^oats a resenblance to sinilar 

ratios appearing In the expression to evaluate f. 

A definition of the resonance utilization, f_, is 

r 



resonance absorption In fuel 



resonance absorption in fuel scattering '’absorption" 

by moderator 

(39) 



The resonance absorption by the uranium is 0^ 
from a thei^al flux absorption analogy i where is an ap~ 
propriately chosen average cross section of uranium over the 
entire resonance flux region. 

Aastiming that the flux over the resonance region is l/E 
dependent, the average cross section is 
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El 

where In ^ * 5*6 for uranium metal (5» p* 273) • The presence 
‘"2 

of the moderator and the geosietric nature of the heterogeneous 

lattice rcquli'es that an effective resonance lnto;;ral CO~ ) 

1*6 

be used for the value ofG^. Equation 4-0 is then 
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where (1, p. 93) 



cr'^ « / ( cr,,) eff ~ = 9.25 + 2l;.7 5/M 

re 



(ij2) 



end s/m la the fuel surface to mass ratio in cm /g. 



finally 
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The scattering ”absorption” by the moderator would be 
^ta from a thornml flux absorption analogy, 



where 






/ p 2^ ^ 



f - 

J X 



s K 
dK 



-37r 



(41^) 



E 



With Fquations 39, 43, and 44 the laagnltudo of the exponent 
of p in Equation 38 is 



V 0 
u ^u 



r ^ 
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(45) 



Murray (1, p. 94) gives an cisiplrioal expresoion for the ef- 
fective resonance inverse diffusion length of uranlura. 



= 0 . 0222 ^^ 



^ is expressed in g/enr^. 



( 46 ) 



where 
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The effective resonance inverse diffusion length for 
graphite is. 



if 



re 



E 




E 

t 



ik7) 



where Z " la the actual tronsport cross section of the graph- 
t 

Ite averaged over the entire resonance flux region* To 
-6 



evaluate the relationship 



T— E - ■ . 

(1 - 00« 9 ) 



(4B) 



is derived from (5» P« 98) 

2 \/'i- - cos e (tt9) 

where cos' '9 « 2/3 A ( 50 ) 



and A is the moderator mss nuraber. 

The method of calculating f is now borrowed to find f^^. 
Equation 19 riodifled is 



V T? 

1 « 1 + P + (E - 1) 

u ^0 



(51) 



Tho values of and were deterrslnoc! from Equations 4^6 

u g 

and 47. Zq and were calculated from Equations and 44 » 
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Using the above values of 4f and F was determined 

^ 8 



from Equation 20 and (E - 1) was determined from Equation 2l|.. 
Table 2 lists all the important lattice material constants. 

Solving Equation ^1 and using the value of f^^ in Equa- 
tion determined the exponent of p in Equation 38 • The res- 
onance escape probabilities for the various configurations. 
Method A, are listed in Table 3* The values of p for the wet 
and dry conf igiirations were the same using Method A. 

2. Method B 

Volkoff and Rumsoy (2, p. 293) express the equation for 
the resonance escape probability as 



where f^ is the resonance utillaatlon of uranium for reso- 
nance neutrons. Noting the form similarity of the exponent 
of Equation $2 to the thermal utilization equations, the fac- 
tor, f„, is 



p axp 



r 

1 - f 



r 



( 52 ) 



f^ « fo (1 + 6) 



( 53 ) 



with 
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0 




( 54 ) 
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changing the cross sections to and Z® to in 

Fquations 31, 33, and 35 gives n solution to Lcuation 
6 is given by Equation 30. 

Substituting tho individual terms in Lquation 54 gives 



^ - 1 » 2 + Y 



(55) 



whore 



V a 









* (“) < ^=1 ''.1 * ’'p * '' w ) 



(56) 



and 



Y= «„V 







_ 1 _ 


(E) 2p Vp + r? Vg 


1 Z’w''w 



^ I 



> (L - 1) 

(57) 



Combining Equations 52, 53, and 55 gives 



p = exp 



. a ■> 6) 

Z + Y - 6 



(58) 



The value of B was again obtained from Equation 24. Tlie 
average cross section of liranium over the resonance region, 
Zq, is given by Equation 43* The scattering "absorption” 
cross section of the moderator over the resonance region, 
is given by Equation 44* Using values from Tables 1 and 2 
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and Equations 20, 2i|., 30, li3, lUl, U6, li.7, 56, 57, and 56 the 
roaonanc© escape probability was deternincd for Method B. 

Those values of p aro listed In Table 3 for tho lattice con- 
figurations considered. The dry values for p were obtained 
by ollalnating tho terms for water in all equations. 

i?’. Buckling 

The theoretical buckling was coraputod using Equations 1 
and 2. The diffusion length for each lattice configuration 
was calculated by (5, p* 200) 

T 2 a ^ « l/h - f) (59) 

Za ® 

where L is tho diffusion length In pur© moderator. The val- 
S 

ue of f used in each case was tho theoretically detertnlncd 
value listed tn Table 3. The Fcrai age, r , of the thermal 
neutrons for tho lattice cojifigtjratlons was assusjed to be the 
eamo aa 7* for the ctoderator. Hurray (1, p. 123) gives the 

p 

value of '>' as 36 ) 4 . cm . The value of kgo for each coaflgiu’a- 
tion was deterciined from tho values of 6 , , f, and p Hated 

in Table 3 for each configuration. 

The volume of a cubical critical reactor havin^j the given 
lattice conflgaration is 
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The square or the lottlco difruniou length end the buckling 
for each of the lattice conf iguj'ations connldered are listed 
in Table 4. volusi® and cube side length of a cubical 

critical reactor for each considered lattice configuration 
are also listed in Table 4. 
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IV . EXPUJ NTA L I?iVF c?T I G A TION 

A. Deaoription of Equipment 

A picture of the subcrltical assembly with a 8.5 in. 
lattice spacing is shown in f'igure 3» The physical dimen- 
sions of the graphite assembly were 60 in. x 60 in. x 79 in. 
Ihe actual graphite used in the construction of the assembly 
was originally 7 in. diameter solid graphite rods. Those 
rods were squared off to a 6 in. x 6 in, cross section (3«5 
in. radius on the corners) and stacked as shown in Figure 3« 
The top five rows of the graphite rods in the assembly wore 
5 in. X 6 in., and the bottom nine rows were 6 in. x 6 in. 

The density of the graphite was determined from a weight 
and dimension measurement of a mass of graphite. The graph- 
ite density was found to be 1.56 g/cm^. 

The graphite assembly was covered on the top and four 
sides with a 10 mil sheet of cadmium sandwiched between a 
0,375 in. thickness of plywood and 0.125 in, of masonite. 

The cadmium sheet gave an effective ”black wall” for thermal 
neutrons. Adding a 0.25 in. total spacing allowance for the 
graphite rods a horisontel section between the ’’black walls'* 
of the assembly was 60,5 in, x 60.5 in. in area. The assem- 
bly was moxmted on a concrete base, and tanks of water were 
placed under the graphite structure to insure that adequate 
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shielding was effected with the neutron source In place. 
Pigxu'e 3 shows tho relative placement of these tanks, 

7he lattice spacing was varied by changing the arrange- 
ment of the uranixim filled procesa tubes in tho assembly. 

The 6 in, lattice, Lattice I, was tho oonfiguratlon with all 
the lattice holes filled vdth procesa tubes. Tho 8.5 in. 
lattice. Lattice II, was t^io configuration as shown on Fig- 
vir© 3» By removing alternate rows of the procesa tubes from 
tho 8.5 in. lattice, the 12 in. lattice. Lattice III, was 
realized. 

Five plutonium-berylllumi neutron sources contained in 
individual right cylindrical containers were placed under the 
center of the graphite assembly. Each of these sources omit- 
ted approximately 1.63 x 10^ neutrons per second j therefore 

6 

the total source strength was approximately 8.l5 x 10 neu- 
trons per second, A cruciform goomotry was chosen for the 
placement of tho five so'urces. One source was at each of the 
cruciform ends, and one was placed at the cruciform center. 
The total source was in this way effectively contained in a 
circle of about 3*5 in. in diameter. 

The cylindrical uranixon slugs were 1.080 in, in diameter 
and 8,40 in. In length. Assuming a IpO 5*111 cladding and a 
200 mil slug ond cap of 2S aluminum, the actual uranium fuel 
size was a rod 1.00 in. in diameter and 6.00 in. in length. 
The density of the uranium was determined from the weights 
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and dimensions of several random sample slu£s. The density 
of natural uranium determined in this way was 19»0 p/cm^. 

The smallest lattice spacing inveotleated, tho 6 in. 
lattice, required ll7 sections of 6061-T6 (cotnncrclal desig- 
nation 61S) aluminum tubing. The tubing oections had an 
outer dianoter of 1.375 in., a wall tlilcknosa of 35 mils, 
and a length of 62 in. Tho comnoaltion of 613 aluminum is 
0.25?^ copper, 0.60J^ silicon, 1.00?^ iruagnealura, 0.255^ chromi- 
um, and tho remainder aluminum. 

The spacing wire composition was 23 aliimlnum and had a 
0*102 in. diameter. The wire was used to center the uranium 
slugs in the process tube. Seven iiranium slugs were loaded 
in each tub© with a helical wrapping of wire around the slugs 
for spacing in the tube. Approximately 10 feet of wire was 
used for each seven slugs. 

Figure l|. shows a scale drawing of an actual unit cell 
of the subci'ltical assembly with a loaded process tube in 
place. ?Tumber seven rubber stoppers were used at both ends 
of the process tubing to contain the water coolant. Tho 
stoppers were also used for the dry configurations to help 
contain the iu*anium slugs during handling of tho loaded proc- 
ess tubes. Allowing two process tube lengths of wire for 
each tub®, the effective thickness of the process tube was 
calculated to be 40 alls. The effective thickness of the 
water annulus was then found to be 0.107 in. 
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Indium foils were mounted on flat aluminum plates which 
in turn were placed on an aluminum holder. The holdei* was 
designed so that the depth of penetration of each foil into 
the slots In the assembly could be accurately detornined. 

The Indium foils were 1«5 in. x 1.0 In. x 0.003 in. and had 
an average weight of 0.59U3 grama. The average thickness 
of these foils was 6.11; ing/ca . Figure 5 shows the foil 
holder and foil placeHSsnt during countln.^. The ^-“oometry of 
the foil with respect to the rlass wall counting tube was 
held constant by use of the count ei' shelf holder shown In 
Figure 5* 

A model 181A scaler manufactured by Nuclear-Chicago and 
a glass wall Gleger tube were used to coxmt the beta activity 
of the irradiated indium foils. A conventional stop watch 
was used for the timing of the counting period. 

B. Procedure 

The uranium filled process tubes were loaded into the 
graphite assembly in varying arrangements for the three geo- 
metric lattices investigated. The annulus between the slug 
and the tube wall was filled water for the wet configura 

tlons. 

The neutron flux distribution was determined by the indi 
um foil activation method. By determining the activity of a 
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foil after a long period of irradiation and making an al- 
lowance for decay during th® process of counting, tho satiara- 
tion activity for a indiun foil was obtained. 

e (61) 

where A,, is the saturation act.ivity, 
io the. activity at timo t, 

^ is the decay constant for indium, 

and t is tho time after removal from the neutron fliix. 

Tho saturation activity is proportional to th© neutron 
flux at an equilibrium, state. The irradiation tin© of tho 
indium foils for each run was in ell cases at least 6.5 hours. 
This length of tin© corres^jonds to more than seven half-lives 
of Indium as the half life for indium Is minutes. The 
indlxim foils used were not exactly all the sano size and 
weight, ?md therefore a norssilising correction was made. Di- 
viding the satiu*ation activity for e given foil by its tfeight 
resulted in a normalized saturation activity. Since a con- 
stant geometry was maintained in tho counter, this operation 
put all of the individual foil activities on tho sane basis. 
This normalized saturation activity is the activity which io 
tabulated In the Appendix, Table 5 through Table 10, for all 
of tho foil positions in each of the lattice configurations 
investigated. 
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Upon romovlnfi a Toll, a throe minute wilnimun decay tiiae 
was given to each foil before beginning the count. This 
period of time elinlnatod any short lived decay components 
in the activated indlusi. All of the counts taken of the 
Indian foils wore three minutes In duration. Tlio decay time 
was taken as the difference between the time of the foil's 
removal from the neutron flux and the nid~tliao of the count 
duration. 

The foil holder slot positions on the front face of the 
subcritical assembly are shown on Figure 6. The horizontal 
X direction spacing was given letter designations A through 
J, The vertical z direction spacing was nxissbered vertically 
1 throiigh 13. The actual distance in in, from the east face 
was the horizontal y direction distance designation. The 
unit of spacing was 6 in. in all cases with the sole excep- 
tion of a vertical spacing change to 5 in. at the position 
whore the lattice structure changed to a 5 in. x 6 in. area. 

A glass wall Glegor tube operating at 9$0 volts was 
used for all of the foil activity deteriuinations. Before 
and after oadi aeries of runs, an operational check on the 
tube and scaler was nade with a sample of radioactive sti'onti 
uia. Tlie background activity remained essentially constant at 
an average of counts per minute during the aeries of co\mt 
ing runs. Dead time corrections were significant only on a 
few of the higher counts. 
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The surv«ya taken of the various lattice configurations 
were as follows: 

Lattices I and II, wet and dry cases, 

a. Vertical at x = 27 in.; y » 18 in. and y » 30 in. 

b. Horizontal, x direction, at y = 18 in. and 
y = 30 in.; z = 30 in. 

c. Horizontal, y direction, at x » 27 in.; 
z ss 30 in. 

Lattice III, wet and dry cesoa. 

a. Vertical at x « 27 in.; y ~ 30 in. 

Iwatticc without uraniuia. 

a. Vertical and horizontal, x direction, at 
y = 30 in. 

b. Horizontal, y direction, at x » 27 in.; 
z = 30 in. 

Figure 6 is a diagram showing the foil positions. 

Glasston© (5;» p. 281;.) shows that the fltix variation in 
the vertical direction at distances not too near the top of 
a rectangular subcrltical assembly shape is 
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The tom Is designated aa the material buckling. Tho 

m 

diraensions of a critical reactor can b© obtained by setting 
P ? 

Bjg * geonetric buckling. V/ith the material buck- 

2 

ling equal to the geometric bucklln.. , B is defined as tho 
buckling. The baae dimensions a and b include the extrapo- 
lation length. At large distances from the assumed source 
plane the harracnic cox^rectiona are assumed to be negligible. 

As the neutron flux is proportional to the indium foil 
saturation activity, a linear plot of th© logarithm of jJ^(z) 
vex'BUS tho ves'tical position z yielded tho quantity - )T 
Yjaich is the slops of this plot. The base dimensions were 
equal In tho assenbly used in this investigation. By adding 
an extrapolation length of 0. ?1 A ^ to each side of the ef- 
fective neutron ”black wall'*, tho square base dimnslons be- 
came 62 in. Equation 63 can now be rewritten as 

= 2 (l] ^ ^ 795 X 10“^ - (6U) 

2 -2 

Where the ixnlts of B are cm , 

Equation 6U was used to evaluate the buckling for th© 
wet end dry cases of all lattice geometries considered. Tho 
slope, - , was measured from the x « 27 in. and y = 30 in. 

positions of the vertical survey plots of each eonf Igxii'ation. 
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V. RESULTS 



The saturation activitloe for given foil positions of the 
coni'iguratlons considered are tabulate<3 lii Tables 6 throtigh 10 
of the Appendix. Plots of these values, figure 7 through Pig- 
uro 16, show the vertical and horizontal distribution of the 
neutron flux in the suhcritical asses’ibly for the lattice con- 
figurations investigated. 

The theoretical and the experisiental values of the buck- 
ling and critical reactor size for the various lattice con- 
figurations considered are Hated in Table 4. Figure 17 shows 
th© relationship between th® theoretical and experlraental 
buckling values for tho wet and <ix*y cases and th© lattice spac- 
ing. 

The y =s 30 in. position of the vertical fl'ux survey plots 
was used for the exporiiaental doteraination of buckling. The 
vertical flux distribution In th© s dii'ection for the y 18 
in. position, Figure 7 and Figure 9, was plotted only to show 
that tho flux variation at this position was also exponential 
in nature. Due to th© critical nature of the slope of these 
vertical flux plots with respxjct to the buckling va3,uc, the 
buckling fro:n the y = 18 in. position was not calculated. The 
buckling, rquation 61}.» is a function of th© square of the 
slope, - , Therefore a very stitall change in the quantity 

- y reflects a much larger chsiige in the value of buckling. 
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A deviation frora a true exponential natiu?e was noted in 
the vertical flux plots at the lower and higher end positions 
of all surveys in the s direction. fho lower end deviation 
was due to excess fast neutrons, 5 to 12 Ifev, fror;! the plu- 
tonima-berylllum neutron source. All of the neutrons from 
the source in this lower portion of the assembly had not as 
yot been moderated, T>5e neutron absorption reaction in indium 
is for neutrons below 1.5 electron volts. The higher ond 
deviation occurred approxlnatoly at the 59 in, height. This 
is the point inhere the lattice spacing chanoee foi* the three 
different geometric lattices considered due to the construc- 
tion of the assessbly. The slope of the flux plots increased 
near this general position in all of the vertical survey 
plots. Statistical error In counting the small amount of 
activity at high a positions accounted for the relative in- 
crease in the dispersion of the counts at tlieso positions. 
These effects can be observed on the plots shown on Figure 7 
through Flgiire 11, 

Ti^e vortical flux for the v;et configuration was greater 
than the vertical flux for the dry conf iguration up to 2 » 

36 in, for Lattice I, Figure 7, at the y = 16 in, position. 

A similar situation existed for Lattice I, Fig\ire 8, at the 
y = 30 in, position with the cross over at a « 50 in. 

Figure 9 shows that the wet-dry cross over point of the 
flux plot for fjattice II at y = 18 In. occurred at z - 30 in. 
Similarly again, the wet-dry cross over point of the vertical 
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flux distribution plot at y * 30 in* for Lattlca II, Figure 
10, and Lattice III, Pisure 11, occurred at a « in. and 
B a 38 in* respectively. The effect of the noderatlon of 
the water was greater than the ex'fect of the water absorp- 
tion for approxisaately the lower half of the aubcritlcal as- 
sejnbly in all lattices investlj-.ated. The water moderation 
effect increased the therml neutron flux by slowing up the 
fast neutrons* In the lower portions of the asaei^ly, these 
fast neutrons are prtmrily from the source. This moderation 
effect by the water in the process tubes >.^s more pronounced 
at the y « 30 In. position than at the y « 18 In. position 
due to the base center position of the source. Over the cen- 
ter, y == 30 in., the neutron flux consisted of raoro fast neu- 
trons at a given low horizontal plane t!ian other positions 
further off center. Tnia excess of fast neutrons was moderat- 
ed in part by the water. At the y « 18 in. position a water 
moderation effect was also experienced by the fast neutrons, 
but tho greater grapliite moderation duo to the greater dis- 
tance from tho source caused the number of fast neutrons to 
be smaller than at tho center positions. The number of fast 
neutrons due to fission was also larger at the center posi- 
tions than at the outer positions of the assembly. TTiis is 
duo to the sine distribution of tho thermal neutron flux in a 
hoz'izontal plane in the assembly. Thermal neutrons caused 
tho fissions in the mt'oral uraniun which produced fast neu- 
trons and other fission products. Tho position of the maxlmus 



thomal noutron flux vms therefore also the position of the 
naxinura fast noutron flux duo to fission. 

The vertical flux survey for the case of the graphite 
lattice without uranivua is shown on P’lgure 32. The dry con- 
figuration flux values for Lattices I, II, and III were cross 
plotted to indicate the cross over point on the flux plot 
without iironiun. The flux plot of T.attice III crossed tho no 
uraniuci lattice plot at a slightly higher point than the flux 
plot of Lattice II did. The fl^ix of Lattice I iiad the lowest 
cross over point of the thro© lattices. The increaeo of neu- 
tron flux due to sx&critical njultiplication in rAttice I ap- 
peared to bo less than the absorption effect of the uraniun. 
This effoct was not nearly so pronounced in Lattices II and 
III. The neutron fltix plot for the no uraniusi lattice ap- 
peared to be trtily exponential in character for all but tho 
low and high extrene positions in the assembly. 

The slopes of tho vertical flux plots for tho x ~ 27 in. 
and y = 30 in. position wore deteralned, and it was found 
that tho nagnitudo of these slopes varied as follows; tho 
flux plot In Lattice II had the snalleat olope and the flux 
plot in Lattice III had the greatest slope. The slope of tho 
flxix plot in Lattice I was slightly greater tiian the Lattice 
III slope. Tho greater the slope of the fliix plot, the siaal- 
ler the valvi© of the buckling becomes. This slope variation 
was true for both the wet and dry com'igurations and was di- 
rectly reflected in tho buckling values calculated. 
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Thoro was little differonco in the effect of subcritlcal 
multiplication between the dry cases of r^attice II and III. 
From a vertical height of hS upwards, the flux in 

dry tiftttlce II was slightly higher than in dry liattice III. 

The absorption difference in these two lattices was apparent- 
ly small. 

The horisontal surveys for symmetry in the x direction 
for z « 30 in. plotted in Figiir® 13 through Figure 15. 

An extrapolation length was added to each end beyond the ef- 
fective neutron ’’black wall”. All of the symetry plots were 
plotted as sine curves. 

The horizontal x direction Lattice I flxix. Figure 13, for 
the wet case wis greater than the flux for the dry ease at 
both y « 18 In. and y « 30 in. positions. Hie horizontal x 
direction r^attice II flux. Figure 13, showed that there was 
very little difference beti/een the wet and the dry cases for 
cither the y » 18 in. or y = 30 in. positions. The y = 30 in. 
position flxix for both lattices was greater in magnitude than 
the y » 16 in. position flux due to the sine distribution In 
the y direction. 

The effect of the water moderation appeared to be greater 
than the water absorption effect in Lattice I at both the 
y « 18 in. and y » 30 in. positions for the x direction survey. 
In the X direction survey for Lattice II, these effects were 
approximately equal for both y positions. 



Lattice II i’liix In the x direction at z « 30 In. had In 
all cases aliaost an ideal sine distribution. Tfae flux in 
wet Lattice I had so'no raridon scatter ovox' the center section 
of the horizontal survey. This variation perhaps was due to 
water noderation having a greeter offoct than tho water ab- 
sorption over the neai' center section. This net effect was 
more pronounced directly over the sotirce. 

Figure 14 shows a horizontal x direction comparison at 
2 * 30 in, of the wet, dry, and without \iraniuai conf insratlons* 
Tjettlco I values of flux for both the wet and dry cases were 
below the flux values in the no uranium Ijattice, The flux in 
wet Lattice I was again higher than the flux in dry Lattice 
I. In Lattice I, the absorption effect at the z » 30 in. 
position was obviously greater than the subcritical jculti- 
pllcaticn effect. Lattice II values of flxix for the wet and 
dry cases arc just slightly below the no tiraniun case. In 
this lattice the subcritical multiplication effect and the ab- 
sorption effects wore approximtely equal, 

A lattice comparison for the horizontal x direction at 
a « 30 in. is shown on Figure 1 $, The Lattice II flux was 
greater than the Ijattice I flux in both the wet and dry con- 
figurations. 

Synnetry (sine distribution) was shown In th® y direc- 
tion on Figure I6 for Iinttlcea I and II at the x « 27 in, and 
z a> 30 in. position. The plot shows that there was little 
difference between the wet and dry oases of either Lattice I 
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or II. In tiattlc© I there was an apparent dirreaslon froa 
ft sine distribution at the center positions of the wet con- 
figuration. Also the slno curve plotted for synaotry vms 
slightly to the left of the points plotted on Plgxire 16. 

This was probably due to the placcnent of the foil holders 
or aoae spreading of the graphite blocks frotrj a uabtice I 
loading. The flux plots of the wet and dry coafigurotlons 
in the y direction for Lattice II are approxlciately the sane 
as the flux plot of the no uronlun lattice. There la again 
aono dispersion at the center position in this comparison. 

The results of the experimental end theoretical bucklings 
determined for all configurations considered are shown on Fig- 
ure 17, The two theoretical methods used to determine the 
buckling for the dry conf igtiratlona showed reasonable agree- 
ment, but the wet case theory was in disagreement. The re- 
sults from Method A were high In comparison with the results 
from Kethod B for wet Lattices I, II and III. 

The experimental buckling values for the wet configura- 
tions wore lovrer than the ex;>erlmental buckling values for 
the dry configurations in all three lattice spaclngs Investi- 
gated. The theoretical values of buckling determined for the 
dry conJ'igurations, ?iethods A and B, wore a little higher 
than the experimental values of buckling for the dry configura- 
tions. Tlio theoretical values of buckling obtained i'rcm 
Kethod A for the wet conf iguratlons were high in comparison 
with the experimental buckling values. The theoretical Method 
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B gave buckling values for the wet conf Igurations wh ch were 
In excollont arreeraent with the exp erlitiental buckling values 
for the wet lattices. 

For the natural uraniua graphite assembly a naxlriuri (!lcl 
occur for the wot and dry cases of th© buckling versus lattice 
spacing curve. Figure 17. l?hls naxisria occui’rcd in both tho 
experimental and theoretical analyses. The C/U aton ratios 
were 73 for Lattice I, ll|L for Lattice II, and 292 for Lattice 
III. Davenport (8, p, 31^) concluded that the .•naxlsiiua buck- 
ling for a given natural uraaluni f-acl assedbly should occur 
at a value of C/U between 50 and 100 * The naximum buckling 
in this Invoatlgation of three lattice spaclxigs occurred at 
about a C/U value of ll^d for both the wot and dry configura- 
tions* 



VI. DI3CU3SI02J OF RESULTS 



Tho theoretical methods employed had the following prob- 
able errors. Accxxrate detailed neutron scattering and capture 
cross section data are not available in the literature for all 
materials. Corrections for varying densities and ptirlties of 
the materials were not confuted, A calculated disadvantage 
factor for uranium was used instead of a measured value. Val- 
ues from the literature were used for tho Fermi age and tho 
various lattice diffusion lengths where experimental values 
again would have given more accurate results for the calcu- 
lated buckling. There were cracks between the graphite blocks 
with a process txibe uranium load in the assembly. These air 
spaces and the air effect around the process tube itself were 
assumed to be negligible. The effect of tesperature and hu- 
midity on the results was also assumed to be negligible. Fi- 
nally, one group neutron theory was used throughout the analy- 
sis. 

The following comments can be made concerning the experi- 
mental methods employed. The harmonic corrections to the flux 
distributions obtained were assumed to be negligible. The 
positioning of the foils and the normalizing of the foil ac- 
tivity were subject to h'oman error. Finally, and most impor- 
tant, the slope of tho In 0{z) versus z direction curve must 
be determined with great accxiraoy. This slop© value is very 
critical in determining the actual buckling. Repeat runs for 
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conslstancy at a given position should bo taken from which 
an average value of activities at a position could be used 

for plottlnr ard the buckling calculation. Varying the foil 

4 

weight at a given position during the runs would tend to can- 
cel out any nornalizing error. 

The feasibility of an actual reactor with the process 
tubing, natural uranium, slug geometry, and water coolant used 
in this investigation was determined. Assuming 1250 channels, 
a water coolant, and 0.^05 in. flow area per channel the mass 
flow is 1.68 X 10^ lb /hr. Assuming a temperature difference 
of 70®F and using Cp =* 1.06 for water at 400®?, the heat rate 
transfer is 366 megawatts. This heat rate transfer is large 
enough for a potential power reactor. The above assumptions 
are conservative in comparison with the experimental critical 
reactor volume determined for the wet 8.5 in. lattice config- 
uration. 

The process tubing did spread the graphite blocks a small 
amount due to the fact that some of the holes were smaller 
than the normal average designed size. A power drill of an 
appropriate size could enlarge these holes so that the process 
tubing would slide in and out of the asserrhly with ease and 
no spread would occur. It was found that a lattice change was 
more easily accomplished using the process tube uranium con- 
tainment than with just the loose slugs. One end of the proc- 
ess tube was beveled to facilitate loading of the wire wrapped 
slugs. It is suggested that a number seven rubber stopper be 
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wall eeatod In the boyeled end and a larger stopper be used 
at the other end. Coolant leakage would in this way be eliral- 
nated* 

The E-3 position foil holder hole was & through hole in 
the y direction. This through hole would be of raore use if 
it were moved up to the E-8 or F-9 positions. The E-3 posi- 
tion was too close to the source for a good syaacnetry flux 
survey. 

The 8.5 in. lattice is the optlmura lattice spacing which 
la possible with the subcritlcal assembly used. Suggested 
ideas for future study are the investigation of other coolants 
and a more detailed analysis of the 8.5 in. lattice for both 
the wet and dry configurations. Also, the effect of the re- 
flection by the water contained below the source inay be in- 
vestigated to give a naxlmum source flux with the available 
sources in plaoo. With the use of an additional imior process 
tube the effect of a variable water annulus could also be in- 
vestigated. 
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VII. coi»cLu:;ioss 

The irwderation effect of the water was groatoi* tiian the 
absorption effoot of the water for approxiHiately the lower 
half of the subcrltlcal assembly for the three lattice spae- 
inga investigated. The thermal neutron flux in the lower 
portion of the assenibly was greater in the wet configurations 
than in the dry cojxfiguratlons, and the revei ae situation 
existed in the upper portion of the asaewbly. The 6 in* wet 
lattice had the lowest overall flux because the absorption 
effect of the greater amount of water and la’aniua was larger 
than the subcrltlcal multiplication effect in this lattice. 

The 12 in. lattice was '‘over moderated” with an apparent small 
subcritioal mltlpllcatlon effect. The optimum lattice in- 
vestigated was the 8.5 in. lattice. In this lattice the water 
had very little effect on the flux distribution due to the 
balance of the water moderation and absorption effects. 

The syraraetry surveys showed that in the x and y direc- 
tions that the flux did follow a sine distribution for both 
the wet and dry conf Iguratlons. 

The two theoretical methods used to determine the buck- 
ling for the dry conf Iguratlons gave resulto with reasonable 
agreement. These theoretical bucklinga for the dry cases were 
high In comparison witii the experimental bucklings determined 
for the dry oases except for the 12 In. lattice. The two 
theoretical methods used to determine the buckling for the 
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wet configurations gave results which were In dleagroonont. 

The method of Murray gave values of the buckling for the wet 
conflgiiratlon of the 6 in, and 8.5 in. lattices which were 
very high in comparison with tho experimental buckling values. 
The method of Rumsey and Volkoff gave values of the buck- 
ling for the wet configuration of tiie 6 in. and 8.5 in. lat- 
tices which were in excellent agreement with the experimental 
buckling values determined. The experimental bucklitig value 
dotornined for the 12 in. wet lattice waa betv;een tho two 
theoretical buckling values determined for this lattice. The 
experimental buckling values for the di*y cases were higher 
than the values for the wet cases in all lattice spaclngs 
investigated. The only experimental positive value of buck- 
ling for tlie wet conf iguratlons was the buckling value for 
the 8.5 in. lattice. 

A maximum in the buckling versus lattice spacing curve 
did occur for both the wet and dry conf igiiratlons. This 
maximum in buckling was found in both the theoretical and the 
experimental values determined. These buckling curves for 
the wet and dry configurations did not cross; hence a '’fall- 
safe" behavior upon the loss of water coolant could not be 
expected from these results. 
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Table 1 



Dlraonsions for a unit cell 
Lattice I (6 In. ) . II (8.5 In.), and III (12 in.) 



* 


uranium rod radius 






1.270 


cm 


"al » 


thickneas of aluminum 


slug 


can 


0.102 


cm 




effective thickneas of 


water annulus 0.273 


cn 


S ' 


effective thickness of 


process tube 


0.102 


cn 


''air 


, effective thickness of 


air 


annulus 


o,kS!^ 


cm 


^1 » 


equivalent inner radius of 


graphite 


2.20 


cm 


Fquival ont 


cross section areas for 


one 


cell 






Uranium 






5.06 


era 


Slug 


can 






0.81; 


cm 


Water 


coolant 






2.61 


cn 


Process tube 






1.06 


-„2 

cm 


Air 








5.65 


ca^ 


Graphite, lattice I 






217.0 




Graphite, I>attice II 






U 34.0 


ca^ 


Graphite, Lattice III 






668.0 


cm^ 


Volume per 


sluf!: 










\ • 


uratiium 






103.0 


ca^ 


Li • 


slug can and cap 






23.0 


cm^ 




water 






55.7 


enr^ 


’'p ' 


process tube 






22.6 


cpP 


'^81 ' 


graphite, Lattice I 






)|625 






graphite. Lattice II 






9250 




’'83 ’ 


graphite. Lattice III 






18,500 


cm^ 
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Table 1 (Continued) 



Lattice I, r^, equivalent outer radius of graphite 8.60 cm 

Lattice II, r 2 » equivalent outer radius of graphite 11. 94 era 

Lattice III, equivalent outer radius of 

graphite 16.76 era 



Table 2 

Lattice cell raeterial constants 



Symbol 


Value 


Source 


^al 


0.01323 


(5) and Eqxiation 25 




0. 01556 cm”^ 


(5) » (11) and Equation 25 


2a" 


0.324 cm**^ 


(1) end Equation 25 


2" 


0.399 cra"^ 


(1) and Equation 25 


2o 


0.0967 cra”^ 


Equation 43 


F 


0.00036 cm”^ 


(5) 


2f 


0.375 cra"^ 


( 1 ) and Equation 25 




0.01057 cra*^ 


(1) and Fqitation 4^i 


2i 


0.354 cra”^ 


(1) and Equation 49 




0.017 cra"^ 


(5) 




0.675 cra“^ 


Equation 26 




0.01992 cm**^ 


(5) and Equation 27 
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Table 2 


( Cont Inued ) 


Sywbol 


Value 


Source 




0,31|72 cra"^ 


(5) and Fquafcion 27 


'u 


O.U218 cm"*^ 


Equation 46 


^8 


0.1059 caT^ 


Equation 47 


cos ^ 


0.0555 


Equation 50 


V'-Jg 


20 


( 3 ) 




0.158 


(5) 


Poriai age, T' 

B 


364 cra^ 


(1) 


X? 

8 


2520 CKi^ 


(1) 


Q 


1.56 g/oh^ 


Measured value 


Cu 


19.0 g/c5!T^ 


Measured value 


S/Vi 


.0828 ca^/gm 


Calculated value. Equation 42 
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Table 3 



Theoretical nuclear lattice conetents 



Lattice 


I 


Method A 


II 


Ill 


DRY 


f 


0.926 




0.870 


0.767 


P 


0.823 




0.909 


0.955 




1.023 




1.062 


0.964 


WET 




0.901 




0.847 


0.751 


P 


0.823 




0.909 


0.955 


* 


0.996 




1.035 


0.963 


Lattice 


I 


Method B 


II 


III 


PRY 


r 


0.920 




0.863 


0.759 


p 


0.836 




0.915 


0.958 




1.034 




1.061 


0.976 




f 


0.882 




0.816 


0.710 


p 


0.632 




0.911 


0.954 


k ^ 


0.965 




1.008 


0.910 


^Evaluated 


losing 6 « 


1.027 and 


1 y * 1.308. 
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Buckling and critical reactor size 







Theoretical 


Method A 




Lattice 




I 


II 


III 


DRY Value 










I,2c.,2 




186 


328 


567 


E ^cr.*2 
m 




41 X 10**^ 


86 X 10 "^ 


-17 I 10*^ 






611 X 10^ 


200 X 10^ 


- 


Cube side, 


ft. 


27.7 


19.2 


- 


WTT Value 










r 2 Z 

L cn 




250 


386 


627 


2 ••2 
la 




-6.5 X 10"^ 


4^6 X 10“^ 


-37 X 10 “^ 


Vycra^ 




«■» 


520 X 10^ 


- 


Cube aide. 


ft. 


- 


26.4 


- 






Theoretical 


Method B 




Lattice 




I 


II 


Ill 


DHY Value 










r2 2 

L csi 




201 


345 


607 


^2 -2 
o ra 




59 X 10**^ 


63 X 10-^ 


-25 X 10”^ 


V,yCW^ 




357 X 10^ 


212 X 10^ 


- 


Cube aide, 


ft. 


23.2 


19.6 


- 


WLT Value 










T 2 2 

L cm 




296 


459 


730 


2 •“2 

cm 
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-23 X 10*^ 


9.7 X 10 "^ 


-85 X 10“^ 


V^cci^ 




MP 


5340 X 10 ^ 


- 


Cube side. 


ft. 


- 


57.4 


- 





Ml 


lit 






*W^46 


Sa»14 


-->jr 


m 


w#« 


%« UA 




• 


%/fl 




41 ,aU« «4M 

«*Wf life 




<•» 


Hi 


V«% 






Si.« 








- 




» » • • «J ^ 
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saflft 
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^C»fi 


•■-V 
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Vfia 






*i» 


M^ 
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Nl 
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V«Ot/k 
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Table Ij. (Continued) 


Lattice 


Experinental 

I II III 


DRY Value 




n 2-2 
b cm 
m 

V.pOa^ 


35 X 10"^ 70 X -5 X 10”^ 

778 X 10^ 275 X 10^ 


Cube aide. 


ft. 30.1 21.3 


WFT Value 




B ^ca”^ 
m 

VjCnP 


-23 X 10~^ 12 X 10"^ -U8 X 10”^ 

3880 X 10^ 


Cuba aide. 


ft. - 51.5 
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A1ISN3Q NOdinaN 



Figuro 1, IfeutrcMi distribution in a uranium grapbito altunlnura lattice 
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Plguro 3. Subcrltic&l aasembly 
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al wire spacer 
0.102 in diameter 




2S al slug can 
0.040 in. 



6 in. 






-6 in. 



61 S al process tube , 
0.035 in. 



air 



water annulus, 
0.112 in. 




3 5 in. radius 




Figure 1}-. Actual unit cell. 
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Figure 7. Vertical flux survey. Lattice I, 
X = 27 in., y = 18 in. 
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Figure 8. Vertical flux survey. Lattice I, 
X s= 27 in., y = 30 in. 
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Figure 9. Vertical flux survey. Lattice II, 

X « 27 in., y =s 18 in. 
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Ftgixro 10. Vertical flux aiurvey. Lattice II, 
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Figiire 11. Vertical flvix survey. Lattice III, 
X s 27 in., y = 30 in. 
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X a 27 in,, y = 30 In. 
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Figure 13. Horl cental flux survey, x direction, Lattice I and 
II, 2 = 30 in., y 18 in. and y » 30 in. 
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Pigur® 15. Horizontal survey, x direction, tiSttlce 



comparison, y « z »30 in. 



SATURATION ACTIVITY c/m SATURATION ACTIVITY c/m 



75 




coraparlacHi, x * 27 In#, z » 30 In. 
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Vertical flux sxirvey lattice I 



Saturation Activity 

0 /n 


Saturation Activity 

oA 


Position 


DRY 


WET 


Position 


DRY 


WET 


El-18 


3903 


4350 


El-30 


6508 


7160 


E2-18 


301-0 


3395 


E2-30 


45 8 0 


5110 


E3-18 


2192 


2460 


E3-30 


3170 


3493 


FU-lS 


1530 


1661 


F4-30 


2085 


2250 


E5-18 


1029 


1095 


E5-30 


1290 


1482 


E6-18 


716 


707 


E6-30 


888 


870 


K7-10 


465 


CO 


F7-30 


552 


620 


E8-18 


335 


313 


K8-30 


376 


378 


K9-18 


238 


230 


E9-30 


268 


253 


r 10-18 


135 


115 


FlO-30 


163 


158 


Bll-lO 


107 


91 


Fll-30 


129 


146 


E12-18 


59 


62 


E12-30 


69 


74 


E13-18 


28 


32 


E13-30 


4i^ 


56 



Table 6 



Horizontal flux survey Lattice I 



Saturation Activity 
C/H 


Sattiration Activity 
C/M 


Position 


PRY 


wifr 


Position 


D?Y 


WET 


a5-i0 


208 


226 


A5-30 


256 


321 


B5-18 


489 


I469 


B5-30 


605 


65 0 


C5~l8 


678 


765 


C5-30 


958 


990 


D!?-18 


893 


950 


P5-30 


1223 


1237 


E5-18 


1029 


1095 


p:5-3o 


1290 


1402 


/5-18 


1030 


1145 


F5-30 


1330 


1453 


G5-10 


952 


1020 


05-30 


1163 


1270 


Hi;-18 


7o5 


748 


K5-30 


939 


994 


1^-18 


505 


504 


15-30 


629 


671 


J5-10 


214 


223 


<15-30 


272 


282 


E5-D 


115 


73 


E5-36 


1306 


ll{,65 


K5-6 


412 


442 


15-42 


1131 


1131 


E5-12 


759 


709 


£5-48 


875 


914 


F5-18 


1029 


1095 


E5-54 


525 


583 


25-214. 


1247 


1320 


E5-60 


180 


150 
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Table 7 



Vertical fl’juc siirvey Tjattice II 



Saturation Activity 

c/m 


t 

a 


Jaturation Activity 

c/n 


Position 


DRY 


WET 


POSitiO‘5 


DRY 


WET 


El-18 




4900 


il-30 


7400 


7815 


E2-18 


366^ 


3002 


E2-30 


5505 


5690 


E3-18 


2720 


2790 


E3-30 


3768 


3975 


EU-18 


1899 


1930 


E4-3o 


2680 


2607 


E5-18 


1241 


1235 


C5-30 


1563 


1571 


E6-18 


801 


784 


E6-30 


1009 


1049 


E7-18 


536 


499 


"7-30 


707 


65l 


E8-18 


376 


332 


E8-30 


4li0 


4S3 


E9~18 


255 


208 


E9-30 


282 


278 


ElO-18 


135 


147 


LI 0-30 


175 


180 


Ell-18 


111 


90 


Ell-30 


135 


llS 


E12-18 


79 


61 


E12-30 


100 


96 


E13-18 


35 


38 


K13-30 


50 


46 
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Table 6 

Horizontal flux swvey lattice II 



Saturation Activity 
C/H 


Saturation Activity 

c/n 


roaition 


DHT 


WET 


Position 


Dir/ 


V/iiT 


A5-1S 


262 


252 


a5-3o 


340 


304 


n5-i8 


632 


607 


B5-30 


791 


728 


C5-18 


090 


926 


05-30 


1170 


1169 


d5-i8 


1105 


1055 


P5-30 


1456 


1435 


E'5-18 


12l}.l 


1235 


E5-30 


1503 


1571 




1230 


1227 


P5-30 


1578 


1550 


G5-18 


1111 


1156 


G5-30 


1487 


1460 


K5-13 


891 


886 


H5-30 


1120 


1189 


15-lS 


602 


595 


15-30 


762 


759 


J5-18 


266 


221 


J5-30 


306 


291 


E5-0 


117 


11^8 


e5-36 


1641 


1572 


E5-6 


511 


541 


E5-42 


1347 


1360 


L5-12 


887 


874 


125-48 


965 


1002 


E5-18 


12U1 


1235 


e5-54 


576 


535 


E5-2U 


1550 


1499 


e5-6o 


195 


160 
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Table 9 

Vortical flux s\irvey Lattice III 



Saturation Activity 
CjM 



Position 


DRY 


WET 


El-30 


7500 


7985 


E2-30 


5945 


6150 


E3-30 


3960 


3990 


EI4.-30 


2813 


2742 


E5-30 


1660 


1658 


E6-3O 


1152 


1093 


E7-30 


669 


637 


EO-30 


450 


425 


F9-30 


265 


234 


Elo-30 


184 


153 


Ell-30 


148 


113 


£12-30 


79 


51 


EI3-3O 


49 


56 
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I^Bttico flxjx survey 


vilthout 


uranium 


Position 


Setiiration Activity 
C/M 


Position 


.Oftturation Activity 
C/H 


El-30 


8365 


1 

0 

1 


172 


E2-30 


6290 


E5-6 


568 


E3-30 


l|li 08 


B’5-12 


992 


Ell-30 


2795 


e5-i8 


1339 


F5-30 


1655 


E5-21| 


1590 


F6-30 


965 


25-36 


1639 


E7-30 


587 


£5-42 


1320 


e3-30 


361 


e5-48 


990 


K9-30 


176 


E5-54 


615 


ElO-30 


128 


e5-6o 


183 


m-30 


81 






£12-30 


54 






1 : 13-30 


22 






A5-30 


320 


G5-30 


I486 


B5-30 


736 


H5-30 


1143 


C5-30 


1161 


15-30 


735 


D5'-3c 


i50D 


J5-30 


318 


’^^^-30 


1669 
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